Ray-wave correspondence in the nonlinear description of stadium-cavity lasers 
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We show that the solution of fully nonlinear lasing equations for stadium cavities exhibits a highly 
directional emission pattern. This directionality can be well explained by a ray-dynamical model, 
where the dominant ray-escape dynamics is governed by the unstable manifolds of the unstable short 
periodic orbits for the stadium cavity. Investigating the cold-cavity modes relevant for the lasing, 
we found that all of the high-Q modes have the emission directionality corresponding to that of the 
ray-dynamical model. 
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Establishing a correspondence between the 
ray/classical picture and the wave/quantum pic- 
ture has been a fundamental problem in the field of 
wave/quantum chaos T]. One encounters this problem 
when trying to understand the emission properties from 
two-dimensional microcavity lasers. In such lasers, as a 
way to extract highly directional emission, it has been 
proposed to deform the cavity shape smoothly from 
perfect circularity 2, 3, 4, 5, 6J. The result is that rays 
start to exhibit a variety of dynamics from integrable to 
strongly chaotic, which is tunable by the deformation. 

The ray picture has been providing a simple and in- 
tuitive method to explain experimental observations of 
emission directionality. For instance, emission direction- 
ality has been associated with the existence of a periodic 
orbit with a particular geometry jQ, J^, drastic shape- 
dependence of emission directionality has been success- 
fully explained by the difference of phase space structure 
, and the far- field intensity patterns have been closely 
reproduced by ray-tracing simulations (|i ^ 0| . 

Among various cavity shapes, the stadium is a sim- 
ple geometry for which ray dynamics has been rigorously 
proven to become strongly chaotic . For almost all ini- 
tial conditions, a ray trajectory explores the entire phase 
space uniformly. Even for such a strongly chaotic cav- 
ity, if one considers refractive emission of light due to the 
dielectric nature of the cavity, the emission pattern can 
become highly directional. Namely, strongly chaotic dy- 
namics and highly directional emission are compatible, 
as was demonstrated by Schwefel et al. Q, who associ- 
ated this property with escape dynamics dominated by 
flow in phase space along the unstable manifolds of the 
unstable short periodic orbits of a chaotic system. 

In this paper, we report further evidence for the ability 
of a ray-dynamical model to describe the lasing states of 
two-dimensional microcavities. Earlier work has focused 
on establishing a relationship between the ray model and 
a few quasi-bound state solutions of the linear wave equa- 
tion, without pumping or gain. Which modes to choose 
for comparison in this case has an intrinsic arbitrariness, 
although plausibility arguments can be made based on 
their Q values. Here we show that the solution of the full 
nonlinear lasing equations for a stadium cavity, uniquely 



determined by the pumping conditions, has highly direc- 
tional emission in good agreement with the ray model. 
This is one of the first pieces of evidence that the multi- 
mode solutions of nonlinear wave equations can be un- 
derstood in terms of the classical limit of its linear coun- 
terpart. 

Firstly, it is shown that a ray model for the stadium 
cavity exhibits highly directional far-field emission. In 
this paper, we fix the aspect ratio of the stadium to 
a/r = 6/7, where r and a are the radius of the semicircles 
and the half-length of the straight segments, respectively 
(Fig. nia))- In the ray model, the appearance of strong 
directionality depends heavily on the value of the refrac- 
tive index of the cavity ni„. We set ni„ = 1.49, which 
corresponds to the index for polymer cavities^. 

The ray model is constructed on the basis of Snell's 
and Fresnel's laws Iilllllili[lS|ll|ll|ll. Inside 
the cavity, the dynamics of a ray is viewed as the motion 
of a point particle moving freely except for reflections at 
the cavity boundary. The ray dynamics can be reduced 
to a two-dimensional area-preserving mapping by intro- 
ducing the Birkhoff coordinates (s,sin0), where s is the 
arc-length along the cavity boundary and cj) is the the 
angle of incidence (Fig. da)). Taking into account the 
dielectric nature of the cavity, we consider the emission 
of rays to the outside of the cavity, which is done in the 
following manner. Each ray is initially assigned a certain 
amount of intensity. This intensity decreases whenever a 
ray collides with the cavity boundary, where the amount 
of the emitted-ray intensity is determined by Fresnel's 
law, while the ray's emission angle is given by Snell's 
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FIG. 1: (a) Geometry of the stadium cavity, (b) the rectan- 
gular unstable periodic orbit. 
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FIG. 2: Ray model simulation: intensity distribution of emit- 
ted rays, I{s,sm(j)). The black solid curves are the unstable 
manifolds of the rectangular unstable periodic orbit marked 
by X. The red curve is a set of points giving the far-field emis- 
sion ai 9 = 210°. The black dashed lines indicate the critical 
lines for total internal reflection deflned by sm(j) = ±l/ni„. 



law. 

In the ray model simulations, we prepare the initial en- 
semble of rays to be uniformly distributed in the phase 
space spanned by the BirkhofF coordinates. After some 
transient, the total intensity of rays inside the cavity de- 
creases exponentially as a function of time . In such a 
stationary regime, we measure the intensity distribution 
for the emitted-rays /(s,sin0). 

The qualitative explanation for strong emission direc- 
tionality in such a chaotic system was given by Schwe- 
fel et al 0; the unstable short periodic orbits act like 
anisotropic "scattering centers" in phase space, causing 
directional flow along their unstable manifolds until the 
critical angle for escape is reached. For shapes like the 
stadium there exists a "line of constant far-field" corre- 
sponding to the set of values of the angle of incidence 
(j) and angular position on the boundary s at which a 
ray refracts in the same far-field angular direction. The 
line of constant far-field closest to the unstable manifold 
then predicts the dominant emission directionality. Note 
that all of the relevant short orbits have closely nested 
unstable manifolds so this gives a unique prediction @. 
The data in Fig 2 are in agreement with this picture. 
In Fig. 2, /(s,sin0) is plotted overlaid with the unsta- 
ble manifolds of the rectangular unstable periodic orbit 
(Fig. 1(b)) located near the critical line for total internal 
reflection, i.e., sine/) — l/riin. For — 1.49, we see that 
high escape intensity regions consist of narrow stripes 
corresponding closely to the unstable manifolds, and that 
two stripes of high intensity (marked by arrows) are al- 
most parallel to the constant far-field curve for 6 = 210° 
(red curve), where the angular coordinate is defined as 
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FIG. 3: The far-field intensity patterns for the stadium cav- 
ity with Uin = 1.49. (A) The ray model (green) vs. the 
Schrodinger-Bloch model (red). (B) The ray model (green) vs. 
the Schrodinger-Bloch model (red: r = 67.77, blue: r = 33.88, 
orange: r = 16.94). 



the counterclockwise angle from the a;-axis. By the four- 
fold symmetry of the stadium each of the other stripes 
is parallel to one of the curves for 9 = 30°, 150° and 
330°, giving a four-fold symmetric directional emission 
pattern peaked around these angular directions. We plot 
the far-field intensity pattern !F{d) in Fig. O which can 
be calculated from the intensity distribution /(s,sin(/)) 
via. J^{9) — J J dsd{sm(j)) I {s, sin (j)) S{tjj{s, (j)) — 9), where 
^{s,(j)) is the emission angle (measured from the x-axis) 
at s for a ray with the angle of incidence (p This is 
the result of the ray model. 

Now, we investigate whether one can find the same 
emission directionality for the lasing states for this cav- 
ity. We describe the light field by the Maxwell equations, 
and assume that the active medium consists of two-level 
atoms obeying the Bloch equations. In the description 
using the Maxwell and Bloch equations, the field vari- 
ables oscillate rapidly at a frequency close to the transi- 
tion frequency loq of the two-level atoms. To perform the 
long-term time-evolutions necessary to obtain stationary 
lasing solutions, we employ the Schrodinger-Bloch (SB) 
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model, which describes the time-evolution of the slowly 
varying envelopes of the field variables 0, 0, 0, Il8| . 
The SB model is given by 



dE 
'dt 
dp 
dt 
dW 
~dt 



-7_lp + kWE, 

~j^^{W -W^)-2k{Ep* + E*p), 



(1) 
(2) 
(3) 



where E{x,y,t) and p{x,y,t) are the slowly varying en- 
velope of the TM electric field and that of the polar- 
ization field, respectively, and W(x,y,t) is the popula- 
tion inversion component. The refractive index n(a;, y) 
is nin{— 1.49) inside the cavity and 1.0 outside it, and 
the linear absorption coefficient a{x,y) is q;l(= const) 
inside the cavity and zero outside it. Space and time 
are made dimensionless by the scale transformations 
(nj„wox/c, n„tJoy/c) {x,y) and Lo^t t. 7^ and 7|| 
are phenomenological relaxation rates, k and p are the 
coupling strength between the light field and the active 
medium, and Woo represents the pumping strength. 

Unless otherwise mentioned, the cavity size is r = 
67.77, for which the perimeter length of the stadium be- 
comes about 105 times as large as the light wavelength 
inside the cavity, i.e., Ai„ — 2-k. The other parame- 
ter values are set as follows: Woo = 0.01, = 10~^, 
7|| = 10~^, ttL = 10~^, n = 0.5, p = 7r/n|^. For the 
above choice of the parameter values, around 100 cavity 
modes have positive linear gain. For such a condition, 
there occur complicated interactions between the modes, 
such as mode-pulling and mode-pushing, that generally 
yield a multi-mode lasing solution [Til w\ . We are inter- 
ested in the far-field patterns arising from these multi- 
mode lasing solutions; we now formulate a convenient 
method for finding this quantity on the basis of the field 
data just outside the cavity. 

Using the Wiener-Khinchin theorem, one can write the 
time-averaged light intensity as 



I(r,t)= / duj Xiui ^|i?J.(r,^)|^ 

1 — >oo ± 
rT/2 



(4) 



where ET{r,uj) = ^ J_^''^^ E{r,t)e"^K In the far-field 

regime (r > 1), we obtain \ET{r,Lu)\'^ — \fT{0,uj)\'^/r, 
where the amplitude fx{0,uj) is given by 



l + i 



'TTh 



<j> dse-''^^^-^"') n,iike + \/)ET{rs,u}), 

(5) 



with k = y^TJn, 



^ 2uj. The integration is performed 
along a closed curve C encircling the cavity, denotes 
a point on C, rig is an unit vector normal to the curve C 
at Ts, and e = r/r = (cos 6*, sin 6*). We define the time- 
averaged far-field pattern !F{9) as the angle-dependent 
part of X(r,i), i.e., 
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FIG. 4: The far-field patterns for the cold-cavity modes with 
even-even parity (solid line) and for the ray model (dashed 
line), (a) C = -0.007054 - 0.002507 i, (b) ^ = 0.009670 - 
0.005197 i, (c) i = -0.006299 - 0.010132 i. Because of the 
pattern's symmetry, only the data in the first quadrant are 
shown. 



In Fig. Ol^A), we plot J-{6) for the stationary lasing 
solution of the SB model. One can see a strikingly good 
agreement between the result from the SB model and 
that from the ray model. It is remarkable that the far- 
field pattern of the SB model reproduces not only the 
highest peaks, but also the tiny ones a±6 5°, 175°, 185° 
and 355°. The magnification of the first quadrant is 
shown in Fig. Ol^B), where the results for the cavity sizes 
r = 33.88 and 16.94 are also presented. Since the increase 
of the r-value results in the decrease of the wavelength, 
one can see that the larger the r-value, the shorter the 
spatial oscillation period becomes. Nevertheless, if we 
average out the oscillations, all these far-field patterns 
show a similar trend with a peak at around 9 = 30°. 
This invariance of the peak position with respect to the 
change of the cavity size convinces us that our system 
with r = 67.77 is well inside the semiclassical regime. 

We note that besides the present study for a/r — 6/7, 
we performed numerical simulations of the SB model for 
various a/r-values ranging from 0.13 to 1.4 while the 
other parameters are set to be the same values as those 
used in this paper. For relatively large a/r-values, we 
could confirm a good ray- wave correspondence. A de- 
tailed study on the condition for the good ray-wave cor- 
respondence will be reported elsewhere. 

One natural approach to explain the appearance of the 
emission directionality in the SB model is based on the 
study of the cold-cavity modes. Analyzing the power 
spectrum for the electric field of the lasing solution for 
r = 67.77, we confirmed that it consists of multiple lasing 
modes, with the major contribution being from six lasing 
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FIG. 5: The dependence of x (far-field emission strength 
around 8 = 30°) on the loss rate — Imf for even-even (cir- 
cle), even-odd (triangle), odd-even (square) and odd-odd (di- 
amond) modes. The modes corresponding to the far-field pat- 
terns in Fig. 13] are marked with arrows. 

modes, and these all exhibit strong far-field emission at 
9 w 30°, 150°, 210° and 330°. From this result, we can 
infer that cold-cavity modes with the above directionality 
are preferentially excited. To investigate the origin of this 
mode selection, we study below the cold-cavity modes, 
focusing on the dependence of the emission directionality 
upon the Q-value, which is one of the important factors 
for the mode selection. 

In the SB model, a cold-cavity mode, E{x, y, t) — 
e~'^^'^ip{x, y) (^ g C), is a solution of Eq. ^ with a = ji = 

0. Namely, V(a;, v) satisfies {v'^ + ^ + 2^ il;(x, y) = 0. 

We plot in Fig. ^a) the far-field pattern of a cold-cavity 
mode having strong far-field emission at « 30°. In the 
gain band of our simulation, i.e., |Re^| < 0.01, we find 
98 cold-cavity modes numerically. As Figs. Efb) and (c) 
show, there also exist cold-cavity modes whose far-field 
patterns are less or even not similar with that of the ray 
model. Such modes, however, turn out to have lower Q 
values as we see below. 

To quantify the far-field emission strength around 

9 = 30° , we compute the quantity x — J20 -^(^) / 



/p X(6') d9, where X(0) is the far-field pattern of the cold- 
cavity mode. Plotting x-values as a function of the loss 
rates — Im^ as shown in Fig. ]E\ we find a clear tendency 
that the lower the loss rate, the stronger the far-field 
emission at 6* « 30°. In other words, at least in this case 
the directional modes based on the unstable manifolds 
also correspond to the high-Q modes and are thus pref- 
erentially selected for lasing, while in general we expect 
that other factors besides the Q value will come into the 
mode selection, such as mode volume and spatial hole- 
burning effects. 

We note that recently Lebental et al. have succeeded 
in experimentally observing highly directional emission 
for a polymer stadium cavity with the aspect ratio a/r = 
0.8 19]. In the experiment, strong far-field emission has 
been observed at 6* « 30°, 150°, 210° and 330°, which 
agrees with our numerical results shown in Fig. |3| al- 
though the aspect ratio for the numerical simulation 
{a/r = 6/7) is slightly different from a/r — 0.8. Car- 
rying out numerical simulations of the SB model also for 
a/r = 0.8, we checked that this slight difference of the 
cavity geometry does not cause a significant change in 
the peak position of the far-field pattern; the peaks are 
shifted by only about 2°. A conspicuous difference be- 
tween the numerical and experimental far-field data is 
that one of the peaks in the experimental data (Fig. 4 
in Ref. ^^) has a three times larger intensity than the 
other three peaks. This is however due to the lifetime 
effect of the laser dye used in the experiment. 

In summary, we demonstrated via numerical simula- 
tions of the nonlinear lasing equations that the stadium- 
cavity laser exhibits a highly directional emission pattern 
in good agreement with the ray model, which predicts 
emission directionality based on the geometry of the un- 
stable manifolds of short periodic orbits. Furthermore we 
confirmed that for the stadium in this parameter range 
all of the high-Q modes exhibit this high emission direc- 
tionality. Further analysis is needed to elucidate to what 
extent this property of the high-Q modes holds when 
one changes the refractive index value, cavity size, cavity 
shape and so on. 
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